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I.  INTRODUCTION 

Systems  ^lre  modelled  in  order  to  understand  and  explain  them  better 
and  53  a prelude  to  action.  Aircraft  dynamics,  for  example,  may  be  identified 
so  tr.at  better  designs  can  be  made,  or  so  that  adaptive  control  actions  can 
be  t5K=r. . Our  attention  in  this  continuing  study  has  been  directed  at 
aspects  of  estimation  and  identification  that  are  connected  with  system 
understanding.  This  study  has  been  aimed  at  two  problem  areas:  multistage 
r.cdeling,  estimation,  a.nd  identification  algorithms  for  dy.namical  systems, 
c:..nd,  state  estimation  and  a parameter  identification  for  a new  class  of 
models  — causal  functional  equations  --  which  describe  wave  propagation  ] 

in  lavered  media  svstems.  i 

i 

The  first  problem  is  concerned  with  developing  estimation  algorithms  | 

both  for  para.meter  and  state  estimation  that  are  recursive  in  the  dimension  j 

of  the  parameter  vector  or  state  vector.  Such  algorithms  will  find  utility  , 

in  system  modeling  work,  where  model  dimension  is  often  a variable. 

The  second  problem  is  concerned  with  develcping  whole  new  theories  of 
state  estimatitr.  and  parameter  ide.ntif ication  for  a new  class  of  equations 
whitr.  we  refer  to  as  causal  functional  equations.  These  equations  are 
co.-.cmucus-time  linear,  ci.me-invariant  wich  multiple  cine  delays.  They 
do  ncc  contain  dari-'atives  or  integrals,  and  no  literature  apparently  exists 
for  c.-.  =m.  Causal  functional  equations  are  applicable  to  diverse  areas  such 
as  reflection  seismclocy,  transmission  lines,  speech  processing,  optical 
thin  coatings  and  EM  problems. 

Curing  the  Summer  of  1977,  Dr.  Keith  Glover  spent  six  weeks  working 
on  this  gra.nt  as  a Visiting  Research  Scholar. 

II.  RESEARCH  PROGRESS 


VJe  have  demonstrated  that  it  is  possible  to  extend  Friedland's 


[A]  bias  estimation  technique,  as  recently  rederived  in  a constructive  manner 


by  Mendel  and  V?ashburn  [1] , to  the  problem  of  estimating  dynamical  states 
and  colored  noise  states  [2] . We  have  shown  how  to  obtain  an  exact  multistage. 
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decomposition  not  only  for  the  state  estimation  equations,  but  also  for 
the  associated  error  covariance  equations.  Additionally,  we  have  obtained 
a second-order  suboptimal  multistage  estimator,  using  a perturbation 
technique.  Whereas  a high-order  matrix  Riccati  equation  must  be  solved  v.’hen 
the  exact  results  are  used,  a matrix  Riccati  equation,  of  the  dimension 
cf  the  colored  noise  states,  must  be  solved  when  th.e  sub-optimal  results 
c.re  used . 

2.  Washburn  [3,  4]  has  generalized  Friedlar.d's  [A]  bias  estimation 
technique  to  partitioned  dynamical  systems.  In  the  general  case,  the 

calculations  are  of  the  dimension  of  the  overall  system,  so  that,  except  { 

for  some  special  but  important  cases,  there  are  no  computational  advantages 
to  the  m.ultistage  approach.  Those  special  cases,  where  there  does  appear 

to  be  computational  adva.ntagas  for  the  multistage  approach  are:  colored  i 

noise  (discussed  above  in  Item  1)  and  weak  coupling  between  the  partitioned  ! 

The  general  results  are  important  in  themselves,  since  they  provide 

J 

the  theory  for  a particular  decomposition  of  the  optimal  state  estimator  | 

for  a system  of  possibly  large  dimension  (i.  c.,  a laroe  scale  system).  ^ 

This  decompositicn  cit’es  added  i.nsight  into  the  structure  and  perforr.ance  of  ■ 

the  minimum,  v.ariance  unbiased  estim.ator.  I.n  addition,  the  methodology  of 
proof  tor  this  multistage  decompositio.n  provides  a mea.ns  for  investigating 
other  decenposi tions  cf  i.nterest. 

Tncer  the  present  grant  we  have  completed  our  study  into  the  development  i 

cf  muitistace  Kalman. 3ucy  filters  for  linear  lumped  parameter  dynamical  systems. 

3.  We  have  developed  time-domain  state  space  models  for  lossless 
layered  media  which  are  described  by  the  wave  equation  and  boundary  conditions 
[5,  i. . Cur  models  are  for  non-equal  one-way  travel  times;  hence,  they  are 
mere  general  than  existing  m.odels  of  layered  media  which  are  usually  for  layers 
of  s-qual  one-way  travel  times.  Full  state  models,  which  involve  2K  states 

for  a K-layer  media  system,  as  well  as  half-state  models,  which  involve  only 
:<  states  have  been  developed  and  related.  Certain  transfer  functions,  which 
appear  in  the  geophysics  literature  in  connection  with  models  of  layered  media 


with  equal  travel  times  have  been  generalized  to  the  situation  of  non-equal 
travel  times.  Our  state  space  models  represent  ii  new  class  of  equations, 
causal  functional  equations,  for  which  we  have  not  Ijcen  able  to  find  any 
literature.  These  equations  are  continuous-time,  linear,  and  contain 
multiple  time-delays.  Their  impulse  response  is  an  infinite  sequence  of 
no.'.-u.-.iiormly  spaced  impulse  functions. 

4.  We  have  prove.n  the  truth  [7,  8]  of  t!ie  following  decomposition  of 
the  solutions  to  the  lossless  wave  equation  in  layer-^c  media:  the  complete 
o_tput  from  a K-layar  media  system,  which  is  comprised  of  the  superposition 
tr  -rimaries,  secondaries,  tertiaries,  etc. , can  be  obtained  from  a single 
state  space  model  of  order  2K  — the  complete  modal  — or  from  an  infinite 
r.u.mber  of  m.od.els,  each  of  order  2K,  the  output  of  the  first  of  which  is  just 
the  primaries,  the  output  of  the  second  of  which  is  just  the  secondaries, 
ett.  This  decomposition  of  the  solution  to  the  lossless  v^^ave  equation  into 


physically  meaningful  constituents  (i.  e.,  prim.aries,  secondaries,  etc.)  is 
called  a canonical  5rem-mer  Series  decomposition,  after  Bremmer,  who  in  1951 
es cablished  a samilar  decomposition  [E] . 

In  many  gecchysacal  situations,  where  reflection  coefficients  are 
qulce  small,  the  decemposition  can  be  truncated  after  secondaries  or  tertiaries; 
he.nca,  it  also  represe.nts  a way  to  approximate  the  solution  to  the  wave  equation  . 
■■'e  have  made  tc.nnecticns  [9j  between  our  state  space  models  and  the 


int-:-ural  equations  given  by  Brammier  [3j  for  generating  the  partial  residuals. 

■•.2  have  shewn  hew  to  go  from  Bremmer 's  integral  equations  to  our  state 
eq'uatitns  by  assuming  ?.  medium  v;ith  a wave  nutdier  that  has  finite  jumps 
(disetntinuities'  which,  occur  at  the  interfaces  and  that  is  constant  within 
a layer.  These  ass'umptions  for  wave  r.u.mber  are  associated  with  what  we  mean 
by  a hcriconcally  layered  homogeneous  earth. 

We  nave  also  demonstrated  that  Bremmer 's  integral  equations  can  be  obtained 
by  the  w.  K.  3.  method  which  gives  approximate  solutions  to  second-order 
differential  equations  [C] . This  justifies  earlier  claims  [7,  8]  that  the 
Bremmer  series  decomposition  can  be  used  to  approximate  the  complete  solution 
to  the  lossless  wave  equation  by  truncating  that  decomposition  after  a small 


number  of  terms. 


5.  Wo  have  developed  [10]  a general  theory  for  describing  reinforced 
events  between  multiple  reflections  in  lossless  layered  media,  which  are 
described  by  tiie  wave  equation  and  boundary  conditions  (e.  g.,  horizontally 
stratified  nonabsorptiva  earth  with  vertically  traveling  plane  compressional 


P.ainforcements  occur  whenever  two  or  more  mv;ltiple  reflections  from 
different  paths  inside  the  media  arrive  at  the  surface  at  the  same  time  so 
that  tnsy  add  (positively  or  negatively)  together.  These  reinforcements 
occur  regardless  of  what  the  travel  tim.e  is  in  each  layer  and  distort  the 
artTr^ranca  of  a seismogram;  for,  they  lead  one  to  believe  that  a significant 
avert  has  occurred  by  the  appearance  of  a large  amplitude  segment  of  the 
seismogram,  whereas,  in  reality,  that  large  event  is  a sum  of  (many) 
smaller  events. 

Our  general  theory  is  applicable  to  a K-layer  m.edia  system  with  non-uniform 
travel  times  and  gives  information  about  the  exact  location  in  time,  number, 
and  amplitude  of  reinforced  events  for  n-aries  (i.  e.,  secondaries,  tertiaries, 
etc.  , where  .n  = 1,  2,  3,  ...  . The  starting  point  for  the  development  of 
this  theory  is  h.endel  ’ s Bre.mmer  series  decomposition  [7,  8]  and  the  operator 
description  of  state  space  m.odels  of  layered  media  [5]  by  means  of  which 
n-ar‘-  reflections  ■' where  n = 1,  2,  3,  ...)  are  generated  and  analyzed 
separaTaly  a.nd  relateu  to  each  other.  Ihe  two  most  significant  multiple 
refl:rCticn3,  saccr.darie? , and  tertiaries,  have  been  studied  extensively.  We 
have  leer nstrated  that  not  only  do  reinforcements  occur  between  the  same  kir\A 
of  multiple  reflections  (o.  g. , between  secondaries),  but  that  reinforcements 
also  ccour  across  different  kinds  of  multiple  reflections  (o.  g.,  between 
secc.-.daries  and  tertiaries)  . 

5.  Because  our  causal  functional  state  space  models  for  a layered 
media  cystem  represe.nt  a new  class  of  equations,  we  have  had  to  study  the 
computer  simulation  of  these  equations.  iivo  computational  methods  have 
been  considered  [11] . In  the  first  approach,  we  discretized  the  time  axis 
and  inserted  states  of  intermediate  delays,  to  arrive  at  a set  of  standard 
finite-difference  equations.  For  our  particular  system,  matrix  multiplications 
can  be  reduced  to  simple  scalar  multiplications.  In  the  second  approach,  we 


defii'.ecl  uiaupiug  rules  for  fVio  tr£insforination  of  states  at  an  interface,  and 
kept  a state  reference  table  for  look-up  and  branching.  The  procedure  is 
similar  to  ray- tracing.  Several  experiments  liave  been  performed  to  sh.ow  the 
trade-off  between  storage  requirement  and  CPU  time-spent  for  the  two  methods. 

7.  V.’e  have  developed  a procedure  for  extracting  reflection  coefficients 
frc".  .‘.oise  data  [12]  which  we  feel  is  a substantial  generalization  of 

sinrlar  procedures  which  have  been  reported  in  the  lirerature  ( [D]  for  example) . 
.-.^scciared  with  these  earlier  procedures  are  Standard  .Assumptions  and  Steps 
•..hith  include  requirements  that  the  data  be  noise  free  and  that  the  observed 
str-sm.ic  data  be  deconvolved.  Our  procedure  avoids  these  restrictive 
recuire.ments . Furthermore,  our  procedure  totally  avoids  the  concepts  of 
z-transforms , minimum,  phase,  spectral  factorization,  forward  and  reverse 
polynomial  manipula cions , etc.,  which  appear  in  the  literature  on  this  subject. 

8.  Several  state  space  realizations  of  seismic  source  signatures  were 
obtained  using  approximate  realization  methods.  This  work  vv’as  performed  by 
nr.  Keith  Glover  who  spent  si.x  weeks  at  the  University  of  Southern  California, 
as  a Visiting  I-.es-arch  Scholar,  during  the  Summer  of  1977. 
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